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Abstract: 

We consider hadronic production of a Z boson in association with a jet and study one- 
and two-loop electroweak logarithmic corrections in the region of high Z-boson transverse 
momentum, ^ Afz; including leading and next-to-leading logarithms. Numerical 
results for the LHC and Tevatron colliders are presented. At the LHC these corrections 
amount to tens of per cent and will be important for interpretation of the measurements. 
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1 Introduction 



Hadronic weak-boson production in association with jets forms one of the most impor- 
tant classes of Standard Model processes. Due to its relatively high cross section and 
clean experimental signature the Z + 1 jet production process can be used for a precise 
determination of the parton distribution functions In particular, the precise measure- 
ment of the Z-boson transverse momentum (px) distribution combined with theoretical 
predictions of correspondingly high accuracy will permit to extract the gluon distribution 
function with unprecendented accuracy at the Large Hadron Collider (LHC). In order to 
achieve desired precision (at the per cent level) radiative corrections should be included 
in the theoretical predictions. 

The next-to-leading order QCD corrections to the total cross section and pT distri- 
bution for hadronic Z + 1 jet production have been known for a relatively long time by 
now [21 IB]- Recent predictions for the Tevatron and the LHC can be found in Refs. 
and respectively. The size of these corrections can amount to several tens of per cent 
depending on the observable under consideration, including jet definition, as well as the 
renormalisation and factorisation scales. 

On the basis of simple arguments regarding the strength of coupling constants one 
might expect electroweak (EW) corrections to be in general much smaller than QCD 
corrections. However, in the region \/I ^ Mw ~ Mz, where the center-of-mass energy 
\/l of partonic scattering becomes much larger than the weak-boson mass scale, the EW 
corrections are strongly enhanced by logarithmic mass singularities. At 0{a^) the leading 
logarithms (LLs) are of the form log^^(s/M^), the next-to-leading logarithms (NLLs) 
are of the form log^^~^(s/M^), etc. These EW logarithms originate from soft /collinear 
emission of virtual EW gauge bosons (7, Z, W) off initial- or final-state particles. Owing 
to the finite weak-boson masses, the real emission of a soft /collinear Z or W boson can be 
observed as a separate process and hence does not need to be included in the definition 
of physical observables. Thus, in contrast to mass singularities in massless gauge theories 
such as QED or QCD, the EW mass singularities of virtual origin are not necessarily 
compensated by corresponding mass singularities from real weak-boson radiation. 

Typically, at \/I ~ 1 TeV the EW logarithms are estimated to yield one-loop correc- 
tions of tens of per cent and two-loop corrections of a few per cent. At one loop it has 
been proven that the EW LLs and NLLs are universal and results applicable to arbitrary 
processes are available E]. For and ZW production at the LHC it was shown 
that these EW logarithmic corrections can become of the order of the QCD corrections 
[H]. Furthermore, in Ref. [HI the EW one-loop corrections to the 7Z production process 
at the LHC were found numerically sizeable, with the logarithmic part providing a major 
contribution to the corrections. The resummation of the EW logarithms was discussed 
first in Ref. [lUj with emphasis on four-fermion processes (for an early discussion see also 
Ref. JJI). General resummation prescriptions for the EW LLs and NLLs were derived in 
Refs. ^11^ and Refs. ^1^^, respectively. For massless four-fermion processes the EW 
logarithms were resummed up to the NNLLs |15j . and it was found that, at the TeV scale, 
the leading and subleading logarithms have similar size and alternating signs, which gives 
rise to large cancellations. These resummation prescriptions were confirmed by various 
two-loop calculations at the LL level 1171 118j and NLL level jl91 l2Uj . Only few two- loop 
results based on a diagrammatic calculation exist beyond the NLL level. In particular, 
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the corrections to a massless fermionic form factor within an abehan massive gauge theory 
j2H I22j were evaluated up to N^LL approximation and agreement with the earlier NNLL 
results was observed. 

In this letter we investigate the Z + 1 jet production at the LHC and the Tevatron. In 
particular, we focus on the high-px region, pT ^ Mz, and discuss the impact of the one- 
and two-loop EW logarithms on the total cross section and the pT distribution for this 
process. We also compare the one-loop logarithmic corrections with the results presented 
in Ref. [22], where the one- loop weak corrections were evaluated numerically. 



2 Analytic results 

The hadronic production of a Z boson with finite transverse momentum px in association 
with a jet receives contributions from the partonic subprocesses 



qq Zg, qg Zq, 
The pt distribution for the hadronic process /11/12 



qg Zq. 
Z + j reads 



dcT^l'^^ 
dpT 



V [ dxi [ dX2 6'(xiX2 - fmm)< 
.3l„;,-^0 Jo 



q=u,d,c,s,b 



fh^,q{xi)fh2,q{x2) + (1 



2) 



dpi 



dpT 



(2) 



where fh,i{x) are the parton distribution functions, fmin = (pt + "T't)^/'S with denoting 



the collider energy and mx = \Jp\ ~\~ M"^. The pT distribution for the unpolarized partonic 
subprocess ij — > Zk reads 



dpT 



PT 



SirNijsli — u\ 



pol 



(3) 



where s = {pi + pj)'^, t = {pi — pz)'^, and u = {pj — pz)'^ are the Mandelstam invariants. 
The momenta Pi.j,k of the partons are assumed to be massless, whereas p^ = for the 
Z-boson momentum. In terms of xi,X2,PT, and the collider energy ^/s, we have 



X1X2S, 



t 



Ml 



[\ — cos ( 



Mi 



u 



-(1 + cos6'). 



(4) 



where cosO = y 1 — Ap^s/{s — corresponds to the cosine of the angle between the 

momenta pi and pz in the partonic center-of-mass frame. The factor 1/Nij in with 



qq — 11^, iiqg — Nc{N^ — 1), and Nc = 3, accounts for the initial-state colour average. 
The unpolarised squared amplitude for the qq Zg process reads 



pol 



1 



tu 



aW + ( ^] aW + ( ^] ^(2) 
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where a and as are the electromagnetic and strong coupUng constants, respectively. The 
amplitude for the g(7-initiated process is easily obtained from @ using crossing symmetry, 



pol pol 



(6) 



The tree-level contribution A^^^ to Q reads 



^^°^=EK)' with /,^. = ^t3^-?^% (7) 



A=L,R 



where T^^ and Yq^ are the weak isospin and hypercharge for left- (A = L) and right-handed 
(A = R) quarks, and we use the shorthands Cw = cos^w and Sw = sin^w for the weak 
mixing angle 

The virtual electroweak corrections were evaluated in the high energy region where all 
invariants are much larger than the weak-boson mass scale, 

|f| > ~ M| for f = s, i, u. (8) 

In this region, using the results of Refs. pi 1191 ITlj . we have computed the one- and two- 
loop electroweak logarithms of the type = log'^dfl/M^) to next-to-leading logarithmic 
accuracy. At L-loop level this approximation includes leading logarithms (LLs) and next- 
to-leading logarithms (NLLs), i.e. contributions of order L|-^ and L?^~^, respectively. 
Contributions of order L?^~^ are systematically omitted. In particular, for the angular- 
dependent logarithms, i.e. logarithms of the type L?^ and L,?^~^ involving invariants f = 
t, u, we use the approximation 

2L _ T 2L , n r l„„ / 1^1 \ T 2L-1 , /n^T 2L"2\ t 2L-1 _ t 2L-1 , /o^T 2L-2\ 



S 



Lf = Lf + 2Llog^ Lf-^ + 0(Lf"^), Lf"^-' = Lf-' + 0{Lf). (9) 



In order to separate the logarithmic electroweak corrections into finite and infrared- 
divergent parts in an elegant way^, it is useful to introduce a fictitious photon mass 
A = M]y. In the present paper we restrict ourselves to the finite part of the corrections 
which is defined as the complete electroweak correction with A = Mw. The remaining 
part, which is not considered here, depends on the gap between the (vanishing) photon 
mass and the weak-boson mass scale and is infrared divergent. This infrared-divergent 
part can be factorized and resummed separately jl2j and needs to be combined with real 
photon radiation. Its size depends on the choice of cuts on the photon energies and angles. 



^ Using the infrared-evolution-equation (TREE) |12j . the electroweak logarithmic corrections can be 
factorized in a SU(2) xU(l)-symmetric and finite part times a U(l)cm-symmetric and infrared-divergent 
part. The SU(2)xU(l)-symmetric part corresponds to the complete electroweak corrections with infrared 
divergences regularized by a fictitious photon mass A = Mw- The U(l)cm-symmetric part contains the 
infrared divergences resulting from the gap between the (vanishing) photon mass and the weak-boson mass 
scale. Within the framework of the IREE these two parts factorize and exponentiate separately. At the 
LL level, the correctness of this IREE-approach has been completely confirmed by process-independent 
two- loop electroweak calculations based on different gauge-fixings |18l 1191 . Recent two- loop calculations 
|19ll2Ull^ have shown that the IREE-approach is applicable also beyond LL level. 
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The small contributions resulting from logarithms of the Z-W mass ratio are neglected. 
Mass-suppressed corrections of order M^/|r| are not considered and, since the contribu- 
tions from longitudinally polarised Z bosons are mass-suppressed in the high energy limit, 
the logarithmic corrections have been computed assuming that the external Z boson is 
transversely polarised. At one-loop level we obtain 



E 

A=L,R 



L| + L? 



Li 



(10) 



where C^^ 



y2 /tA^2 



for quarks, with Cq^ = 



+ Cq^l are the eigenvalues of the electroweak Casimir operator 
3/4 and Cq^ = 0. The above result takes into account LLs and 
NLLs of soft /collinear origin as well as NLLs resulting from the running of the coupling 
constants.^ The latter NLLs are, however, cancelled by the collinear logarithms that are 
associated with the final-state Z boson [Hj and thus do not appear in the result ()1U() . 
In the one-loop approximation and for the process under consideration the electroweak 
corrections can be split into electromagnetic and weak parts in a gauge-invariant way. 
The weak corrections (which result from virtual Z- and T^-boson exchange) are easily 



obtained from through the substitution C[ 
electromagnetic charge of the quark q. 
At two-loop level we obtain 



ew 

q\ 



Qq, where Qq denotes the 
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(11) 



where h = -41/(6c|,) and 62 = 19/(6s|,) are the one-loop /5-function coefficients associ- 
ated with the U(l) and SU(2) couplings, respectively. We note that at the two-loop level 
the purely weak corrections cannot be isolated from the complete electroweak corrections 
in a gauge-invariant way. The LLs as well as the angular-dependent subset of the NLLs in 
i.e. all contributions of the form L| with f = s, t, u, have been derived from Ref. |19j . 
There, by means of a diagrammatic two-loop calculation in the spontaneously broken elec- 
troweak theory, it was shown that such two-loop terms result from the exponentiation of 
the corresponding one- loop corrections. The additional NLLs of the form L| in have 
been obtained via a fixed order expansion of the general (process independent) resummed 
expression proposed in Ref. This NLL resummation Jl] relies on the assumption 

that effects from spontaneous breaking of the SU(2)xU(l) symmetry can be neglected in 
the high energy limit. 

^We note that in addition to the logarithmic corrections in the partonic subprocesses involving b 
quarks receive also NLL corrections proportional to Mt /M^, which originate from the Yukawa interaction 
in the heavy-quark doublet At partonic level these corrections can be relatively large. However, their 
contribution to the hadronic cross section is almost negligible owing to the smallness of the &-quark density 
of protons and anti-protons. We have checked that the correction resulting from these Yukawa terms does 
not exceed one per mille at the LHC. 
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3 Numerical results and discussion 



In this section we discuss numerical effects of the one- and two-loop logarithmic EW 
corrections on the Z + 1 jet cross sections at the LHC and the Tevatron. Results presented 
here are obtained for the following values of the electroweak parameters: a = a{M'^) = 
1/127.9, = 0.231, Mz = 91.19 GeV and Mw = Mz^Jl - s^. When discussing the 
(relative) statistical error, defined as Afistat/c = l/^v/iV with N = £xBR(Z — > /, ui) xcjlO) 
we include all leptonic decays of Z with the branching ratio BR(Z = 30.6% and 

assume total integrated luminosity of £ = 300 fb~^ and C = 11 fb^^ for the LHC [2^ and 
the Tevatron (23) respectively. 

To calculate the hadronic cross sections we use the leading order parton distribution 
functions MRST2001LO at the factorisation scale fMp = p\- Similarly, the strong 
coupling constant as is evaluated (for five active flavours) at the scale /i^ = p\ and, 
consistently with the parton distribution analysis, we use the one- loop expression with 
as(M|) = 0.13 [ini- 

No rapidity cuts are imposed for the purpose of this study, as their effect on the 
distributions does not appear to be relevant, both at the LHC and the Tevatron. 
We checked that at the LHC a rapidity cut \yz\ < 4 decreases the px distribution at 
Pt = 100 GeV only by less than 0.5% and for px ^ 100 GeV the effect of the rapidity cut 
becomes negligible. 

The lowest order (LO) transverse momentum distribution of the outgoing Z boson at 
the LHC is shown in Fig.^. We also display seperately the lowest order contribution from 
the longitudinally polarised Z boson. Its smallness justifies to treat the radiative correc- 
tions to the unpolarised process as if the Z boson would be only transversally polarised. 
The next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) results, i.e. 
one- and two-loop corrected distributions, are also presented. 

The size of the corrections relative to the leading order p^ distribution is shown in 
Fig. ^p. These relative corrections increase with px- For high px they yield a negative 
contribution of tens of per cent to the px distribution in the one-loop case and a positive 
contribution of several per cent in the two-loop case. More precisely, for 1 TeV < px ^ 
2TeV, the logarithmic one-loop electroweak correction decreases the cross section by 24% 
to 39%, and the two- loop correction moves the result up by 3% to 8%. The biggest 
(positive) contribution to the two-loop result is given by the LL terms, with NLL terms 
accounting for about 50% of the size of the LL contribution at high px but with an 
opposite sign. This observation demonstrates that the control of NLL effects in the two- 
loop approximation is required to predict the cross section at a level better than 10% for 
high Px- We note that for the one- loop result the EM contribution grows with px but 
even at the highest px considered it does not exceed 3-4% of the total electroweak one-loop 
logarithmic correction and thus is negligible. 

The size of the one- and two-loop electroweak logarithmic corrections to the px distri- 
bution clearly makes them relevant for future measurements at the LHC. The significance 
of the corrections is also illustrated in Fig. [21 where we plot the relative one- and two- loop 
corrections to the total cross section obtained by integration over px with px > Px^* as a 
function of p™*. Additionally, in the same figure we show the statistical error evaluated 
assuming the LHC luminosity and the branching ratio quoted above. The one-loop cor- 
rection is evidently bigger than the statistical error in the entire range of px considered 
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Figure 1: (a) Transverse momentum distribution for pp Zj at ^/s = 14TeV: LO 
(dotted), NLO (dashed) and NNLO (solid) result for unpolarised Z and LO contribution 
from longitudinally polarised Z (dash-dotted), (b) Relative electroweak correction to the 
lowest order unpolarised pT distribution for pp Zj at ^/s = 14TeV: 1-loop LLs+NLLs 
(dashed), 2- loop LLs (dash-dotted) and 2-loop LLs-|-NLLs (solid). 

here. The two-loop correction is larger than (or comparable to) the anticipated statis- 
tical error for values of pT up to about 1500 GeV. In fact, the correction from the LL 
term is significantly larger, being, however, strongly reduced by the NLL correction. This 
demonstrates that at the two-loop level at least leading- and next-to-leading logarithms 
will be relevant for a reliable prediction. Of course, a full estimate of the experimental 
error requires including systematic effects what leads to an increase in the error value. 
Yet, Fig. [21 indicates relevance of the electroweak corrections for the precise analysis of 
experimental data. 
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Figure 2: Relative electroweak correction and statistical error for the unpolarised in- 
tegrated cross section for pp Zj at ^/s = 14TeV as a function of p™*: (l+2)-loop 
LL+NLL (solid), 2-loop LL (dash-dotted) and 2-loop LL+NLL (dashed) correction and 
statistical error (shaded region) with respect to the lowest order cross section. 

We also compared the lowest order pT distribution and one-loop logarithmic correction 
for the LHC with the results of Ref. >23,;- In their calculations the authors of Ref. [231 
introduced rapidity cuts, chose the scale of running as to be n'^ = s, and while evaluating 
the one-loop corrections did not consider the electromagnetic part but only the weak 
contributions. Taking into account these differences we find a good agreement for the 
lowest order distribution and one-loop corrections, at the level of a per cent for the 
latter. This shows that the main part of the one- loop EW correction is given by the 
logarithmic contribution. 

Results for the Z + 1 jet production at the Tevatron are shown in Fig. |2K,b and in 
Fig.EJfor the pT distribution and the integrated cross section, respectively. Contrary to the 
LHC case, at the Tevatron the one- and two-loop logarithmic electroweak corrections to 
the cross section do not bear much significance for the precise measurement, as illustrated 
in Fig. m Comparing our results with those of Ref. jSB] we note a disagreement for the 
lowest order distribution by roughly a factor of ten at high p^ values. 

4 Conclusions 

Experiments at the LHC will for the first time explore energies in the range well beyond 
1 TeV. In this region, where s ^ ~ M'^, the electroweak radiative corrections be- 
come important due to large logarithms of ,§/M^. In the present paper we have studied 
their impact on hadronic Z-boson production at large transverse momentum. We have 
calculated the one- and two-loop electroweak corrections in logarithmic approximation, in- 
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Figure 3: (a) Transverse momentum distribution for pp Zj at ^/s = 2 TeV: LO 
(dotted), NLO (dashed) and NNLO (solid) result for unpolarised Z and LO contribution 
from longitudinally polarised Z (dash-dotted), (b) Relative electroweak correction to the 
unpolarised lowest order distribution for pp — > Zj at ^/s = 2 TeV: 1-loop LLs+NLLs 
(dashed), 2-loop LLs (dash-dotted) and 2-loop LLs-|-NLLs (solid). 



eluding leading and next-to-leading logarithmic terms. These corrections increase rapidly 
with pt- At the LHC, for values of beyond 1 TeV the one-loop terms amount to several 
tens of per cent, reaching up to 40 per cent at 2 TeV. These results are in good agreement 
with those based on an explicit one-loop calculation [22]. To fully control the prediction 
with a precision that may be reached in the final stage of LHC operation, even the dom- 
inant two-loop terms must be included. Significant compensations between leading and 
next-to-leading two-loop logarithms are observed. The combined two-loop effect of several 
per cent is comparable to the anticipated statistical error at the LHC. Contrary to the 
LHC, at the Tevatron the one- and two-loop logarithmic electroweak corrections do not 
have numerical significance. 
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Figure 4: Relative electroweak correction and statistical error for the unpolarised cross 
section for pp — > Zj at ^/s = 2TeV as a function of p™*: (l+2)-loop LL+NLL (solid), 
2-loop LL (dash-dotted) and 2-loop LL+NLL (dashed) correction and statistical error 
(shaded region) with respect to the lowest order cross section. 
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